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I ■ Abstract. When a thin sheet is crushed into a small three-dimensional volume, it 

'^ ■ invariably forms a structure with a low volume fraction but high resistance to 

Q . further compression. Being a far-from-equilibrium process, forced crumpling is not 

, ^, ! necessarily amenable to a statistical description in which the parameters of the 
initially flat sheet and the final confinement fully specify the resulting crumpled 

►^ state. Instead, the internal geometry and mechanical properties of the crumpled 

^ ■ ball may reflect the history of its preparation. Our X-ray microtomography 



experiments reveal that the internal 3-dimensional geometry of a crumpled ball is 
IT) '■ in many respects isotropic and homogeneous. In these respects, crumpling 

fr^ I recapitulates other classic nonequilibrium problems such as turbulence, where a 

O ■ system driven by long-wavelength, low-symmetry, forcing shows only rather subtle 

fingerprints of the forcing mechanism. However, we find local nematic ordering of 
the sheet into parallel stacks. The layering proceeds radially inwards from the outer 
yj , surface. The extent of this layering increases with the volume fraction, or degree of 

'^ I compression. 



While the global 3-dimensional arrangement of a crumpled object is very complex, 
a remarkable feature of the geometry is that a large fraction of the surface area is 
relatively fiat. The curvature imposed on an elastic sheet by external confinement 
or forces dwells is concentrated largely in a network of ridges [H [2] which meet at 
vertices known as developable cones [3l H]. For a perfectly elastic sheet, the work 
done in crumpling is stored in the elastic energy of these focused deformations, 
which is partitioned in finite fractions of bending and stretching energies [1]. For 
most familiar examples of crumpled sheets such as plastic, paper, or metal foils, the 
strains at ridges exceed the yield strain and the ridges become irreversibly creased 
into folds. In the crumpled regime, these focused structures must interact simply 
due to geometric confinement. Understanding the mechanics of an interacting set of 
folds is a formidable challenge, and it is crucial to obtain experimental insights into 
their 3-dimensional arrangement. 

Some important lessons are learned by studying a lower-dimensional version of 
the problem, that is, the packing of a 1-dimensional elastic curve confined in 2- 
dimensions. Experiments on the packing of a wire [5] in a circular cage, and of 
sheets pulled through a circular hole [B] showed a new feature: as a curve is confined, 
it starts organizing into parallel arcs. Based on a recent lattice simulation of this 
problem [Tj, it has been argued that true long-range nematic order can be obtained 
in the thermodynamic limit. In this lower dimensional problem, it is possible for 
the curve to pack without any stretching, and therefore, no focused structures are 
produced even in the limit of zero thickness. However, 2-dimensional sheets confined 
in a spherical space are subject to greater geometric frustration - they cannot ac- 
commodate the constraints presented by confinement in all three directions solely by 
bending. They are forced to stretch and develop gaussian curvature by bending in 
two directions [8J . It is thus important to come to terms with the fully 3-dimensional 
problem as there are qualitative differences between these situations. 

Two recent sets of numerical simulations of purely elastic [S] and elasto-plastic 
[TU] self-avoiding sheets in 3-dimensions have studied the distribution of elastic energy 
in ridges and the scaling of compressive forces with confinement. The configurations 
shown in Ref. [10] clearly indicate that layering also occurs in 3D. While there is no 
direct spatial quantification of this phenomenon, they find that the energy becomes 
dominated by self-contacts as the confinement increases. 

It is a challenge to experimentally probe the interior geometry of a 3D crumpled 
object. Ref. pi] studied by laser profilometry the statistics of folds and vertices in 
an unfolded crumpled sheet. However, unfolding the sheet leads to a loss of spatial 
information about the interactions of the folds and of the sheet's final crumpled con- 
figuration. A recent paper [T2] studied by laser sheet illumination the early stages 



of energy focusing in an elastic sheet. X-ray tomography is a tool that is well-suited 
for providing structural data in the highly crumpled regime. In recent X-ray tomog- 
raphy experiments [131 E] a 2D analysis of the geometry was performed, in terms of 
tangent vectors to the curves obtained in 2-dimensional slices of the crumpled object. 
This analysis showed parallel orientation of tangent vectors; thus both simulations 
and experiment indicate that layering of facets is a geometric element to be consid- 
ered alongside focus! ed structures such as ridges and vertices. However, a truly 
3-dimensional experimental characterization of the highly-crumpled state remains to 
be done. 

In this article we present measurements by X-ray microtomography of the location 
of a crumpled sheet in a 3D volume, along with fully 3-dimensional analyses of the 
curvature, orientation and layering of the sheet. We concentrate on an analysis the 
radial dependence of these metrics of the geometry, in order to identify any possible 
anisotropy and inhomogeneity imposed on the sheet by the method of crumpling. 

In our experiments, circular aluminum sheets with thickness t = 25.4/im and 
initial radius, Rt are hand-crumpled into approximately spherical balls with final 
radius Ro- We have studied balls with average volume fractions of ~ 6%, 8.5% 
and 22%. X-ray computerized tomography (CT) is used to image the interior of 
the crumpled sheet. A CT scanner (Skyscan 1172) images the radiation from a 
divergent x-ray illumination from a point source. Transmission images are taken 
after 0.3 degrees rotations of the sample over a total of 180 degrees. These images 
are reconstructed into a three-dimensional stack of gray scale images with a typical 
size of ~ 2000 x 2000 x 2000 pixels. The voxels are uniformly distributed in the 
volume, and are cubes with a linear dimension of 8.85/im ?a t/3, thus ensuring that 
the thickness of the sheet is fully resolved. As can be seen in the 2D slices through 
the crumpled balls shown in Fig. [H we are now in a position to extract various 
quantitative measures of the disposition of the sheet in the crumpled volume. 




Figure 1 : Reconstructed slices through an equatorial plane of three crumpled spheres with 
average volume fractions oi (p = 6% {Ri ^ 3.5cm, Ro ~ 0.72cm), 8.5% Ri f« 5cm, Ro ~ 
0.82cm) and 22% {Ri Ki 7cm, Rq ~ 0.75cm), respectively. To remove background noise, all 
reconstructed images are thresholded before analysis. 



1 Results and Analysis 



1.1 Mass Distribution 

The simplest representation of the geometry of the sheet is the radial dependence of 
the volume fraction 0(r). We ensure by calibrating with a sample of known geometry 
that the mass of film determined from the reconstructed images is not affected by 
x-ray absorption through the volume of the sphere (see SI). As the individual blue 
curves in Fig. [2]A. show for nine different crumpled sheets with the same average 
volume fraction (p = 8.5%, there is considerable variation of 0(r) from sample to 
sample. The black curve, which is the average over these samples, more clearly shows 
a trend of volume fraction increasing from the interior to the exterior of the sphere. 
Thus the most elementary analysis reveals a signature of the low-symmetry route 
to the final crumpled state, where the confining forces are radial and inwards. The 
large variability within any of the individual samples reflects the very heterogenous 
distribution of void space with the ball. This heterogeneous distribution has often 
been quantified in terms of a fractal distribution of mass [TU [151 [13 [IS] • However, 
given that there is an overall radial density gradient, a fractal dimension computed 
from an average over the volume of the ball is not a useful measure of heterogeneity. 
The data shown in Fig. [2l3 compare the radial gradient of volume fraction for 
three different degrees of confinement. The functional dependence on radial distance 



does not change significantly with volume fraction over this range of (j). This must of 
course change at extremely low and high [l3j confinement where a more homogeneous 
distribution might be expected. 
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Figure 2: Mass Distribution. A (left), Radial dependence of volume fraction, (j){r) versus 
r/Ro, the radial distance r from the center of mass normalized by Rq, the radius of the hall. 
The center of the sphere is defined by its center of mass. The data shown are for spheres 
with a radially averaged volume fraction of cp = 8.5%. The volume fractions are normalized 
by (j), where the over-bar represents a further average over samples. The blue curves are for 
nine separate samples, to give an indication of the sample-to- sample variability. The black 
curve is the average of these samples. Error bars indicate the standard deviation. The 
average mass distribution increases k. 30% from, the centre of the sphere to the exterior. 
B (right), Normalized volume fraction <p{r)/cj), versus the radial distance from the center 
of mass for spheres with (p = 6%, 8.5%, and 22%. The red open circles and green open 
squares corres! pond to spheres with (p = 6% and 22%o respectively. The black curve 
is the average of 9 spheres with tp = 8.5%. The solid line is a best-fit to the equation 
cpfit = II 3(7/5 (1 + Cr'^ /Rq), with an adjustable parameter C. 



1.2 Orientation 

We next move from the location of the sheet, to the orientation of sheet within the 
volume. From the gray-scale 3D image, we find the surface normal n at all points 
along the surface. The spatial orientation of the surface normal is quantified by the 
direction cosine n ■ f, where f is the unit vector from the centre of the sphere to 
the surface point. The distribution of this direction cosine, shown in Fig. |3JA. is 



nearly uniform, thus indicating a near-isotropic distribution of sheet normals. This 
is a surprising result, given that the crumpling process might be expected to break 
symmetry between radial and azimuthal directions, and perhaps favour orientation of 
the sheets parallel to the confining surface, in onion-like fashion. This expectation of 
alignment by the outer surface is not borne out even as one approaches the surface. 
In Fig. [3|3, we plot the average n ■ f as a function of normalized radial position 
r / Ro- The orientation, h ■ f, does not deviate from the value expected for random 
orientation, except perhaps within a small region close to R = Ro- Except for this 
boundary layer, the orientation of surface normals is both isotropic and homogeneous 
within the volume, despite the sheet's route to the crumpled state. 
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Figure 3: Orientation. We quantify local orientation of the sheet by the dot product h ■ r 
of the surface normal and the radial unit vector from the centre of the hall, as shown in 
the 3D patch of sheet on the left. We display orientation data from spheres with <f) = 6%, 
8.5%, and 22% (red open circles, black diamonds, and green open squares respectively) . The 
8.5% data are an average of 2 spheres. A (left) , Histogram of the sheet orientation, \h- r\. 
All three volume fractions have a nearly isotropic orientation with only a slight tendency 
for alignment of the normals with the radial vector. B (right), Average orientation versus 
distance from the center of mass for (p = 6%, 8.5%, and 22%. The solid line displays the 
average of an isotropic distribution of orientations. 



1.3 Curvature 

The next higher order in the description of the local geometry of the sheet is the 
curvature of the sheet. In principle, the curvature can be determined from a complete 



knowledge of the vector field of surface normals h ^18j . However, numerically finding 
the gradient of this field introduces undesired inaccuracy. We chose to make an 
independent measurement of the curvature by identifying connected patches of the 
sheet, and fitting ellipsoids to these patches. From the fit parameters, we were 
able to obtain the two principal radii of curvature Ri and -R2(> Ri) at every point 
on the patch. The unusual geometry of the crumpled state presents a technical 
challenge here. While there are many regions of low curvature, the curvature can 
change sharply at the stress-condensed regions. Thus the surface has to be fit with 
rather small patches. It was not possible to infer large radii of curvature with great 
numerical accuracy by fitting to small patches in space. We therefore focus our 
analysis on highly-curved regions with radii of curvature < 50t, which can be reliably 
fit. In Fig. HB, we show for three volume fractions, the high-curvature part of the 
histogram of the principal radii of curvature. There is a peak in the histogram at Ri 
~ lot, independent of volume fraction. The radius associated with the peak in the 
histogram, Ri = lOt, and with the largest fitted radii Ri^2 = 50t, are indicated by 
open circles in Fig. HJA.. The value of curvature radius at the peak is much higher than 
that required to introduce plastic folds in the aluminum. The yield stress for Al is 
145MPa [12], and a radius of curvature for yielding of 250t may be inferred from the 
corresponding yield strain. It thus appears that there is a limiting value of curvature 
beyond which the crumpling process creates new features, rather than compressing 
existing folds to sharper dihedral angles. As we discuss below, though, the peak at 
small radius of curvature does not imply a high level of energy condensation. (Also 
see SI for further details). 

From the curvatures we have determined, we identify surface points with one 
and with both radii of curvature below a cutoff radius, and associate these with 
folds and vertices, respectively. The cutoff radius, Tc, for each (p is chosen to be the 
radius at which the histogram falls to half the value of the peak in the curvature 
distribution, as shown by the dotted lines in Fig. HJB. The number of high-curvature 
points identified by applying an arbitrary threshold obviously depends on the value 
of Tc- With the definition we adopt, we find that the fraction of surface points 
associated with folds is 37% and with vertices is 4% for the 0=8.5% sample (Fig. 
I1P,D). The large fraction of surface points with small radii of curvature indicates 
that in this regime of confinement, stress condensation is incomplete and that it is 
not correct to think of the sheet as largely fiat. In Fig. HP we display the fraction 
of points with Ri < re (folds and vertices), as a function of radial position. Fig. UP 
shows the same quantity for points with Ri < R2 < r^ (vertices). It is evident from 
these data that the high-curvature parts of the sheets are homogeneously distributed 
through the sphere, and show no evidence of being preferentially generated at the 



confining walls. A plausible explanation for the radial increase in volume fraction 
might have been that the exterior part of the sphere allows more gentle curvature, 
and therefore that it is energetically preferred to have more of the sheet reside in the 
exterior. In light of the observation of a homogeneous distribution of high-curvature 
features, and of the isotropic orientation of the sheet normals, this explanation is 
not tenable. This homogeneous curvature distribution contrasts with crumpling in 
2-dimensions [6], where curvature is different in the bulk and at the boundary, where 
the object assumes the curvature of the confining wall. Possibly this is due to the fact 
that the spherical boundary condition imposes a gaussian curvature on a fiat sheet, 
which unlike in the 2D analogue, is not possible to accommodate without stretching 
deformations. 
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Figure 4: Curvature. Displayed are measurements of the high curvature regions of spheres 
with (t)=6, 8.5% and 22% (red open circles, black diamonds, and green open squares re- 
spectively). The 8.5% data are an average of 2 spheres. A (top left). A reconstructed 
image from a crumpled hall with (j)=8.5% is shown to give a sense of scale for the radii 
of curvature. B (bottom left). Normalized histogram of the larger and smaller principal 
curvature-radii Ri and R2 in units of t. Radii of curvature between 3t and 50t are shown. 
There is a peak in the histogram at Ri/t ~ 10 for all three 4>. Locations with radius of 
curvature below a cutoff radius, r^ are designated as high- curvature regions; r^ is chosen to 
he the radius at which the normalized histogram falls to half its peak value and its value is 
indicated for each volume fraction by the dashed vertical line. From left to right, the three 
open circles show the radius at the peak in Ri (Wt), re for 8.5% (2At), and the upper limit 
of curvature measurements (50t). C (top right). The fraction of s! urface points with at 
least one high- curvature direction (folds and vertices) versus distance from center of mass. 
D (bottom right) The fraction of surface points with two high- curvature directions (vertices 
only) versus distance from center of mass. The high- curvature regions are homogeneously 
distributed. 



1.4 Nematic Ordering 

Thus far, we have concentrated entirely on local descriptors of the geometry. We now 
turn to the stacking of facets that is evident in Fig. [1] As a quantitative measure 
of this local nematic ordering, we study the correlation of the 3D surface-normals 
along the sheet. From each surface point we search along the normal direction for 
other surfaces in parallel alignment. We label surfaces stacked together in this way 
within a chosen search radius, and count the number of sheets layered together in 
each stack, m. In Fig. [5j\, we show as a function of radial position, r/Ro, the 
fraction of surface points which are participating in stacks composed of m =3 to 
7. These data for 0=8.5% reveal a tendency to stack that increases approximately 
linearly in the radial direction. In Fig. [5B we show the probability of occurrence 
of 3-stacks and 6-stacks for different volume fractions. At all three volume fractions 
we study, the probability of 3-stacks is similar, and increases approximately linearly 
in the radial direction. However, 6-stacks are just beginning to form in the outer 
layers of the sphere at lower volume fractions, and develop more strongly at higher 
volume fractions. The implication is that stacking is initiated from the outside and 
progresses inwards as crumpling proceeds. Thicker stacks could then be produced 
by accretion or folding of thinner stacks, and are not due to a strong increase in the 
overall fraction of layered sheet. 
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Figure 5: Stacking. Stacks are identified by starting at all surface points on the sheet, 
searching in the direction of the normal h, and identifying the number, m, of other surfaces 
with anti-parallel normals within a radius of Kit. For example, m=l represents an isolated 
part of the sheet, and m=4 represents 4 closely-stacked surfaces. The solid circle in Fig. [7] 
gives a visual sense of the search radius. A (left). Fraction of surface points participating 
in stacking versus distance from the centre of mass for different values of m in one sphere 
with (j) = 8.5%. The symbols represent stacks containing m=3,4,5, 6 and 7 facets (darkest 
to lightest colored diamonds respectively) . The stacking fraction increases with r/Ro for all 
values of m. Inset, histogram of surface orientation, \h ■ f\, for points on stacked surfaces 
(i.e. m >2) only (blue crosses). We also show for comparison the histogram for all surface 
points (red plusses). This shows that the stacked facets are also oriented isotropically! . In 
the panels on the right we show stacking fraction versus r/Ro for (j) = 6%, 8. 5%, and 22% 
(red open circles, black diamonds, and green open squares, respectively) for (B) m=3 and 
(C) m=6. As (j) increases, stacks become thicker. 

As mentioned earlier, the crumpled state spontaneously develops structural rigid- 
ity at very low volume fractions without externally imposed design. This rigidity has 
been attributed to the formation of ridges with high buckling strengths[9|fT0]. While 
there is no direct mechanical verification of either mechanism, the layering shown in 
Fig. [5] may also be a contributing factor to mechanical rigidity. That is, the structure 
may be stabilized against external compression by forming multi-layered walls rather 
than pillars. The force threshold to buckle a ridge into smaller ridges [20], or for Eu- 
ler buckling of a planar region, scales with the bending modulus of the sheet, which 
grows as t^, the cube of its thickness. If the sheets in an m-stack do not slide relative 
to each other either due to frictional contact, or by the ridges that might delimit the 
edges of the layered region, then the rigidity of such a stack scales as {mt)^. Stack- 
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ing can thus greatly enhance mechanical rigidity. The relationship between stacking 
and stress condensation remains to be clarified. It has been suggested in studies of 
phyllotaxy in cabbages that layered arrangements are nucleated between stiff leaf 
stems [2^. A similar mechanism may be operative here, with the folds playing the 
role of the stem, and forcing sheets in near-contact to make stacks. Even though the 
stacks are formed under radial compression they are also isotropically arranged, and 
potentially strengthen the structure against forces in any direction. This is shown 
in the inset to Fig. [5JA. where we separately display the distribution of orientation 
with and without the stacked layers. (A different result was obtained when a 2D 
determination of orientation was performed [13]). Any arbitrary force applied to the 
crumpled ball will result in both compression and shear internally. The response of 
the stack to such a generalized stress state may involve both the self-avoidance of 
the sheet as well as friction between the sheets. 

The detailed geometrical characterization of the crumpled state raises questions 
that we are currently attempting to address by simultaneous studies of structure 
and dynamics. However, the principal inference to be drawn from the structural 
information presented here is that several aspects of the geometry are to a good 
approximation homogeneous and isotropic, and therefore perhaps amenable to sta- 
tistical treatment. Despite the low-symmetry path to the crumpled state, it is rather 
remarkable that from the vantage point of a location in the interior of the ball, no 
local measurement of geometry points the way to the exterior of the ball. 

1.5 Methods and materials 

To find the radial distribution of mass, the reconstructed image stack is first thresh- 
olded, and then despeckled to remove detector noise, (^(r) is determined by counting 
occupied voxels in spherical shells of thickness Ar = 3.3t. We ensure by measure- 
ments on a calibration sample that absorption effects do not affect the measured 0(r) 
(see SI). We choose to define the radial distance r from the centre of mass of the 
ball. The exterior hull of the ball is not perfectly spherical, nor is the centre of mass 
necessarily at the geometric centre of this hull. However, at these average volume 
fractions, the differences between these measures is not significant for the results 
we discuss. The final radius Ro is determined from the steep fall-off in 0(r) at the 
exterior of the ball, and is chosen to be where 0(r) k. 0.0450. Ro is consistent with 
measurement of the caliper radius of the crumpled ball, averaged over orientation. 

The orientation of the sheet is determined by applying a 3D canny edge-detector 
to find the surface normal ft at all points. The edges are smoothed by a Gaussian 
filter. We use a simple threshold to fill in low-intensity surface points. 
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To measure curvature, the sheet is first thinned by identifying the number of 
neighbors in a 3 x 3 volume surrounding that point. This procedure reduces both 
surfaces of the sheet to a one-voxel-thick edge. All detected surface points are only 
counted once. The entire sheet surface is decomposed into connected patches that are 
fitted to algebraic functions (see SI) from which the local curvatures are determined. 
Our analysis is optimized toward accurate determinations of the high- curvature fea- 
tures of the sheet. 

We identify the number of facets stacked between two parallel surfaces by search- 
ing in the direction of the normal at each surface point, ni, at a radius of up to ~ 16t, 
(see the solid circle below Fig. HJA. for an visual indication of the search radius). We 
identify all other surface points along this search direction with normals 122 that are 
antiparallel to ni as defined by the criterion rii ■ n2 < —0.97. We then return to the 
raw image to count all occupied voxels between these points with antiparallel nor- 
mals, m-stacks are identified by dividing this count by the sheet thickness measured 
in pixels. When the sheets are very tightly stacked, the interior surfaces in contact 
are not resolved in the images. However, our identification of the number of sheets, 
m, in the stack is not affected by this complication, since we count the entire mass 
between the two outermost sheets in the stack. 

We gratefully acknowledge Whitey Hagadorn and Diane Kelley for access and 
guidance to the microtomography apparatus. We also thank H-q. Wang, G.R. 
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